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CONCENTRATION  QUENCHING  IN  Ndj^Y^_jjPjOi4  CRYSTALS 

John  M.  Flaherty*  and  Richard  C.  Powell 
Department  of  Physics,  Oklahoma  State  University 
Stillwater,  Oklahoma  74074 

ABSTRACT 

An  extensive  Investigation  of  the  spectroscopic  properties  of 
NdjjYj^_^P50j^4  Is  reported.  Time-resolved  site-selection  spectroscopy 
and  photoacoustic  spectroscopy  techniques  were  used  as  well  as  standard 
optical  spectroscopy  measurements  which  were  made  as  a function  of 
Nd  concentration,  temperature,  excitation  wavelength  and  laser  power. 
Results  show  that  Nd  Ions  occupy  a variety  of  slightly  perturbed  sites. 

No  significant  energy  transfer  appears  to  be  taking  place  between  Ions  in 
these  inequivalent  sites.  Additional  experimental  data  presented  here 
suggests  a new  quenching  model  for  neodymium  pentaphosphate  crystals  based 
on  spatial  but  non-spectral  energy  migration  to  surface  quenching  sites. 
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I.  INTRODUCTION 

The  development  of  the  technology  of  fiber  optics  has  resulted  In  a re- 
quirement for  miniature  optically  pumped  lasers  with  low  threshold  and  high 
galn.^  Stoichiometric  laser  materials  such  as  neodymium  pentaphosphate 
(NdP50i4)  meet  these  requirements. Despite  the  significant  amount  of 
work  that  has  been  done  recently  there  are  still  some  aspects  of  the  lumi- 
nescence properties  of  concentrated  Nd  materials  which  are  not  well  under- 
stood. Complete  characterization  of  the  optical  properties  of  this  class 
of  materials  Is  Important  not  only  for  technological  design  considerations 
but  also  for  obtaining  a basic  understanding  of  the  physical  Interaction 
mechanisms  Involved. 

The  most  Important  property  of  stoichiometric  laser  materials  Is  that 
they  exhibit  unusually  weak  concentration  quenching.  This  allows  the  con- 
centration of  neodymium  Ions  In  NdP502^4  to  be  very  high  without  significant 
degradation  of  the  fluorescence  and  lasing  properties  of  the  material.  The 
characteristics  of  the  concentration  quenching  In  stoichiometric  laser  mate- 
rials are  different  from  those  In  other  well  characterized  materials  such  as 
Y^Al50j^2 physical  mechanism  for  the  quenching  process  is  not 
understood  at  this  time.  Developing  a model  to  explain  concentration  quenching 
In  NdP^Oj^^  Is  an  important  problem  In  understanding  this  class  of  materials. 

' To  date  three  basic  mechanisms  have  been  proposed  to  explain  the  con- 
centration quenching  characteristics  of  stoichiometric  laser  materials. 

All  have  been  proven  unsatisfactory.  The  first  mechanism  postulated  was 
cross- relaxation  between  pairs  of  Nd  Ions  In  the  bulk  of  the  crystal  which 
predicted  a quadratic  concentration  dependence  for  the  quenching  rateS-8. 

This  prediction  is  Inconsistent  with  experimental  results. A second  suggestion 
for  a quenching  mechanism  Involved  crystal  field  overlap  mixing  effects  between 
Nd^  lons^^  but  this  has  been  proved  wrong  by  the  observation  that  the 
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oscillator  strengths  of  the  transitions  are  Independent  of  concentration.^^ 

Finally  energy  migration  to  ions  In  randomly  distributed  sites  throughout  the 
bulk  of  the  crystal  acting  as  "sinks"  where  the  energy  Is  dissipated  radlatlon- 
lessly  has  been  proposed. Unfortunately  this  model  requires  the  unreal- 
istic assumption  that  the  concentration  of  Ions  In  quenching  sites  is  the  same  for 
all  samples  Independent  of  total  Nd  concentration.^  We  recently  reported 
preliminary  results  indicating  that  spectral  energy  transfer  between  Nd3+  ions 

in  slightly  different  types  of  sites  is  not  taking  place  in  NdP50i4  as  it 

♦ 

does  in  materials  such  as  Y3AI5O12 even  at  low  concentrations.^^ 

In  this  paper  we  report  the  results  of  an  investigation  of  the  mixed 
crystal  system  of  where  we  have  utilized  some  new  spectroscopic 

techniques  in  an  attempt  to  better  understand  the  concentration  quenching 
mechanism  in  these  materials.  Laser  time-resolved  site-selection  spectroscopy 
was  used  to  investigate  energy  transfer  and  the  results  show  that  migration 
among  Nd3+  ions  in  non-equivalent  sites  does  not  take  place.  Photoacoustic 
spectroscopy  (PAS)  measurements  were  made  to  characterize  radiationless  relax- 
ation processes  among  the  excited  states  of  Nd^*^  in  these  materials.  The 
results  are  not  consistent  with  those  expected  for  PAS  signals  originating 
mainly  from  those  phonons  due  to  relaxation  processes  from  higher  absorption 
levels  to  the  metastable  level.  A model  based  on  surface  quenching  of  the 
excitation  is  proposed  which  can  explain  the  observed  PAS  results  and  is 
consistent  with  all  the  other  spectroscopic  observations  obtained  on  this 
material.  The  possible  mechanisms  for  the  quenching  process  and  the  effects 
of  radiative  reabsorption  are  also  discussed. 

II.  SAMPLES  AND  EXPERIMENIAL  EQUIPMENT 

The  samples  were  grown  from  solution  at  Philips  Laboratory.  Six  samples 
of  IM.Yi.^PsO  2^4  ranging  from  x*1.00  to  x*0.10  were  investigated.  These  con- 
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slated  of  small  crystallites  varying  In  sice  from  about  2.0  mm  on  a side  to 
about  O.SO  na. 

For  optical  spectroscopy  measurements  the  samples  were  mounted  In  a 
cryogenic  refrigerator  to  control  the  temperature.  For  photoacoustic  spectro- 
scoDv  the  samoles  were  mounted  in  an  acoustical  cavity  and  the  signal  detected  by 
« l-in.  electret  microphone.  The  excitation  source  for  obtaining  PAS  and 
fluorescence  excitation  spectra  was  a 1000-W  tungsten/halogen  lamp  passed 
through  a 1/4-m  monochromator.  For  measuring  fluorescence  llfetlsws  and 
obtaining  time-resolved  spectroscopy  data  the  excitation  source  was  a nitrogen 
laser-pumped  tunable  dye  laser  with  rhodamlne  6-G  dye.  This  provided  a pulse 
of  about  10  ns  in  duration  and  less  than  0.5  ^ in  halfwidth.  Sample  emission 
was  analyzed  by  a one-meter  monochromator  and  detected  by  a cooled  RCA 
C31034  photomultiplier  tube.  To  obtain  steady  state  spectral  data  the  Incident 
light  was  chopped  and  the  signal  processed  by  a PAR  lock-in  amplifier.  For 
processing  pulsed  signals  a TAR  boxcar  averager  triggered  by  the  laser  pulse  was 
used.  The  time  gate  of  the  boxcar  was  either  set  to  observe  the  fluorescence 
spectrum  at  a specific  time  after  the  laser  pulse  or  scanned  to  obtain  the  fluor- 
escence as  a function  of  time  at  a specific  wavelength. 

Figure  1 shows  block  diagrams  of  the  experimental  setups  for  time-resolved 

spectroscopy  and  photoacoustic  spectroscopy  measurements. 

III.  GENERAL  SPECTROSCOPY  AND  CONCENTRATION  QUENCHING  RESULTS 

The  absorption  spectrum  of  NdP50j^^  has  been  presented  elsewhere.^ 

It  is  similar  to  the  spectrum  of  Nd^*^  Iona  in  other  hosts  with  lower  doping 

e 

concentrations.  The  largest  absorption  band  occurs  at  approximately  8000  A 
corresponding  to  the  transition  from  the  ^19^2  state  to  a band  consisting 

of  the  unresolved  components  of  the  ^H9/2  ^^5/2  l*vel8  (see  the  energy 

diagram  of  Nd^"^  shown  in  Fig.  2).  The  second  largest  absorption  band  occurs 
near  5800  X and  is  associated  with  the  unresolved  levels  of  the  ^^^^2 
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^Cj/2  l«v«ls.  After  absorption  Into  any  of  the  higher  energy  levels, 
radiationless  relaxation  occurs  to  the  ^^3/2  metastable  state.  Fluorescence 
emission  then  takes  place  from  the  two  crystal  field  components  of 
this  level  to  the  various  Stark  components  of  the  multlplets.  Branching 
ratios  determined  previously^  indicate  that  most  of  the  energy  la  emitted  in 
the  transitions  terminating  on  the  ^lxi/2  ^^9/2  centered  around 

1.06  pm  and  0.89pm,  respectively.  The  fluorescence  spectrum  is  discussed  In 
more  detail  in  the  next  section. 

Figure  3 shows  the  fluorescence  excitation  spectra  of  the  ^^3^2 
monitored  at  8880  X at  room  temperature  for  both  NdP50}^^  and  a mixed  crystal 
of  Ndo. 5^0. 5^5^16*  spectra  are  quite  similar  for  the  pure  and  mixed 

samples.  In  particular,  the  excitation  bands  centered  at  ' 8000  X and  ' 5800  X 
are  comparable.  The  relative  strength  in  the  absorption  spectrum  of  these  two 
bands^  correlated  with  the  fact  that  the  8000  X excitation  is  approximately 
the  s.ime  or  less  intense  than  the  5800  X band  In  the  excitation  spectrum  sug- 
gests that  the  longer  wavelength  band  Is  less  efficient  In  contributing  to  ^^3/2 
fluorescence  than  the  shorter  wavelength  band. 

Figure  4 shows  the  fluorescence  lifetime  of  the  ^^3/2  level  at  room  temper- 
ature as  a function  of  Nd  concentration.  The  lifetime  decreases  from  about 
220  PS  for  the  NdQ^I^Y^  9P5O24  sample  to  about  95  ps  for  NdP50^^.  A concentration 
quenching  rate  can  be  defined  as 

Wq  - - ^o"^  (1) 

where  Tq  Is  the  Nd  fluorescence  lifetime  in  the  absence  of  any  quenching  inter- 
action. This  quenching  rate  is  also  plotted  in  Fig.  4 as  a function  of  con- 
centration. For  an  intrinsic  lifetime  of  * 220  ps  the  plot  of  the  quenching 

annears  to  varv  approximately  linearly  with  concentration. 

The  fluorescence  lifetime  of  the  ^f'3/2  level  was  measured  with  different 
wavelengths  of  excitation  as  a function  of  temperature  between  14  K and  300  K. 
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The  fact  that  a significant  variation  in  lifetime  with  excitation  wavelength 
was  observed  la  discusaed  in  the  next  aectlon.  The  temperature  dependence  of 
the  fluorescence  lifetime  for  NdP^Oj^^  for  two  different  wavelengths  of  exci- 
tation is  shown  in  Fig.  5.  There  is  a small  Increase  in  lifetime  between  14 
and  300  K.  A similar  trend  is  observed  for  other  samples  and  other  excitation 
wavelengths. 

IV.  TIME-RESOLVED  SITE-SELECTION  SPECTROSCOPY  RESULTS 

The  optical  spectral  lines  of  impurity  ions  in  solids  all  exhibit  in- 
homogeneous broadening  to  some  degree.  This  is  due  to  the  fact  that  imperfec- 
tions and  internal  strains  give  rise  to  non-equivalent  crystal  field  sites  for 
the  impurity  ions.  In  some  cases  the  site  differences  are  great  enough  that 
transitions  from  ions  in  different  types  of  sites  can  be  distinctly  resolved. 
High  resolution  lasers  have  been  used  to  selectively  excite  impurity  ions  in 
specific  crystal  field  sites  in  glass  and  crystal  hosts.  By  tuning  the  laser 
excitation  through  an  inhoroogeneously  broadened  absorption  band  it  has  been 
possible  to  observe  significant  variations  in  fluorescence  spectral  parameters 
such  as  peak  positions,  intensities,  and  radiative  and  radiationless  transition 
probabilities.^^  If  pulsed  lasers  are  used,  time-resolved  spectroscopy  tech- 
niques allow  the  observation  of  the  fluorescence  spectrum  at  specific  times 
after  the  laser  pulse.  After  selective  excitation  of  ions  in  specific  sites 
the  fluorescence  spectrum  in  some  systems  has  been  observed  to  evolve  with 
time  into  spectra  exhibited  by  ions  in  other  types  of  sites.  This  implies 
that  energy  transfer  is  occurring  between  ions  in  different  types  of  sites 
and  the  observed  time  dependence  can  be  used  to  characterize  the  properties 
of  the  energy  transfer  and  help  in  identifying  the  mechanism  responsible  for 
the  transfer  process. 

For  investigating  site  selection  spectroscopy  in  the  I^x^l-x^S^lA  ^Y*^*"* 
we  pumped  into  the  absorption  band  consisting  of  the  unresolved  ^^7/2 
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7 ^G^/2  level*  because  this  transition  is  known  to  be  "hypersensitive"  to 

;i  changes  in  the  local  environment  of  the  ion  11,15,16  and  for  other  systems 

such  as  Nd  in  mixed  garnet  crystals  we  found  that  pumping  into  this  level 

resulted  in  easily  observable  site-selection  and  energy  transfer. The 

I 

' fluorescence  was  monitored  from  the  ^^3^2  level  to  the  various  components 

ot  the  ^1^^2  **ti*fe  ma  n told  in  order  to  minimlre  homegeneous  broadening 

it  ‘he  transitions  due  to  rad lat lonless  relaxation  processes.  As  the  narrow 


line  laser  excitation  was  scanned  across  the  broad  absorption  band  at  low 
temperatures,  for  all  of  the  samples  Investigated,  changes  in  the  relative 
fluorescence  peak  intensities,  peak  positions,  and  lifetimes  were  observed. 

1 • ' 

For  example,  for  the  NdPjOj^^  sample  at  14  K the  lifetimes  ranged  from  63us 
to  105  ys  depending  on  the  excitation  wavelength  whereas  at  300  K they  varied 

'i  '■ 

f between  about  102  and  115  vis.  Figure  6 shows  an  example  of  the  different 

A 

fluorescence  spectra  which  are  obtained  for  the  NdP5034  sample  at  14  K for 
|n  pumping  at  different  wavelengths.  The  total  integrated  intensity  of  the 

If  , 

, fluorescence  is  approximately  the  same  for  each  excitation  wavelength  but  the 
[ distribution  of  the  intensity  among  the  five  lines  and  the  fluorescence  llfe- 

I times  differ  significantly.  The  lifetime  is  67  ys  for  5784.5  A excitation  and 

O 

77  ys  for  5837  A excitation.  The  highest  energy  line  (labeled  1 in  the  figure) 


shows  the  greatest  variation  with  pumping  wavelength  but  there  are  relative 
changes  in  the  intensities  of  all  of  the  other  lines  also.  Figure  7 shows  an 
expanded  picture  of  line  1 as  a function  of  excitation  wavelength  at  14  K for 
the  Ndo.5^0.5*’5°  24  crystal.  Here  the  changes  in  line  shape  and  peak  position 
are  quite  evident.  For  the  5831  A excitation  the  line  is  resolved  into  two 
distinct  peaks. 

These  spectral  as  well  as  lifetime  changes  as  a function  of  excitation 
wavelength  can  be  attributed  to  the  selective  excitation  of  Nd  ions  in  non- 
equivalent crystal  field  sites.  If  rapid  energy  migration  occurs  among  the 
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Nd  Ions  the  spectra  should  evolve  with  time  into  some  characteristic  spectrum 
representing  the  average  distribution  of  the  energy  among  all  of  the  Ions  In 
different  types  of  sites.  We  attempted  to  observe  this  by  monitoring  the 
spectra  as  a function  of  time  after  the  laser  pulse  In  the  range  from  1 to  200 
us.  No  significant  time  evolution  of  the  spectra  was  detected  In  this  range 
for  either  of  the  two  samples  investigated,  NdP^Oj^  and  NdQ  ^Yq  5^5014^  for 
any  wavelength  of  excitation  at  any  temperature. 

In  order  to  quantitatively  describe  any  tine  dependent  relative  changes 
in  the  spectra  obtained  for  two  different  excitation  wavelengths  a fluorescence 
deviation  function  can  be  defined.  The  spectrum  obtained  for  a given  excitation 
wavelength  is  represented  by  a five  component  function  each  component  cor- 
responding to  the  ratio  of  the  Integrated  intensity  of  one  of  the  lines  of  the 
transitions  to  the  total  Integrated  intensity  of  all  of  the  lines. 

Thus  each  fluorescence  spectrum  can  be  characterized  by  a five  component 
vector.  Any  correlation  between  the  fluorescence  spectra  for  different  excitation 
wavelengths  can  be  seen  by  forming  the  vector  difference  between  the  functions 
describing  the  spectra  for  each  excitation  wavelength.  Figvire  8 shows  typical 
results  for  the  fluorescence  deviation  as  a function  of  time  after  the  laser 
pulse  for  both  the  NdPjOj^^  sample  and  the  Ndo.5YQ  5^5^14  various  temperatures. 

A value  of  D equal  to  0 would  indicate  an  identical  distribution  of  relative 
intensities  in  the  five  spectral  lines  for  the  two  excitation  wavelengths. 

The  presence  of  energy  transfer  between  Nd  ions  in  different  types  of  sites 
should  result  in  a time  dependent  decrease  in  D toward  zero.  Instead  for  all 
cases  investigated  D remains  within  experimental  error  constant  with  time 
Indicating  that  no  energy  transfer  is  taking  place  between  ions  in  non-equivalent 
sites. 

The  temperature  dependence  of  the  fluorescence  deviation  function  for 
NdPjOj^  is  shown  in  Fig.  9.  The  error  bars  indicate  the  spread  in  the  values 


of  D for  measurements  at  different  times  after  the  laser  pulse  at  each  tem- 
perature. Since  the  decrease  in  D with  temperature  occurs  without  any  associated 
time  dependence,  this  temperature  dependence  can  be  attributed  to  the  reduced 
ability  of  the  laser  tr  selectively  excite  Nd  ions  in  specific  sites  at  high 
temperatures.  This  may  be  a result  of  homogeneous  broadening  of  the  absorp- 
tion levels  due  to  electron-phonon  interactions. 

The  site-selection  spectroscopy  results  presented  above  show  that  Nd  ions 
in  pentaphosphate  hosts  occupy  different  types  of  crystal  field  sites.  Although 
the  exact  nature  of  these  different  sites  is  not  known  they  are  probably  due 
to  differences  in  local  strains  resulting  from  structural  Imperfections  and 
defects  in  the  crystal.  The  similar  relative  intensities  observed  for  different 
excitation  wavelengths  imply  similar  concentrations  for  the  different  types  of 
sites.  The  fact  that  no  energy  transfer  could  be  detected  between  ions  in 
different  types  of  sites  Is  quite  surprising  in  light  of  the  fact  that  similar 
studies  made  on  Nd  ions  in  mixed  garnet  crystals  show  strong  energy  transfer 
between  ions  in  non-equivalent  sites. ^7  it  is  important  to  note  however  that 
spatial  energy  transfer  between  Nd  ions  in  equivalent  sites  may  take  place 
without  spectral  energy  transfer  to  ions  in  different  types  of  sites.  This  type 
of  transfer  would  not  be  observable  by  the  time-resolved  site-selection  spectro- 
scopy techniques  used  here  and  is  discussed  further  in  Sections  6 and  7. 

V.  PHOTOACOUSTIC  SPECTROSCOPY  RESULTS 

Photoacoustic  spectroscopy  techniques  have  recently  been  developed  to 
characterize  radiationless  decay  processes  of  ions  in  crystals. By 
comparing  PAS  results  with  fluorescence  excitation  spectra  it  has  been  possible 
to  determine  dominant  relaxation  channels. 20,21  lqw  resolution  photoacoustic 
spectra  were  obtained  at  room  temperature  on  the  concentrated  neodymium 
pentaphosphate  samples.  The  1000  W tungston-halogen  lamp  chopped  at  110  Hz 
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was  used  as  an  excitation  source.  Figure  10  shows  the  results  obtained  on  the 
NdPjOj^  sample  and  the  NdQ  ^Yq  . The  signal  to  noise  ratio  on  the  more 

lightly  doped  samples  was  too  poor  to  give  meaningful  results.  Note  that  the 
relative  peak  heights  of  the  major  absorption  bands  In  the  PAS  spectrum  are 

comparable  to  those  In  the  absorption  and  excitation  spectra. 

These  PAS  results  are  more  difficult  to  Interpret  than  those  In  the  previous 
samples  investlgated.^®"^^  The  photoacoustic  signal  normalized  for  the  energy 

of  the  Incident  light  beam  at  phase  angle  0 can  be  expressed  as  the  sum  over 
all  the  relaxation  transitions  which  generate  heat^® 

where  Ng  Is  the  number  of  photons  absorbed,  the  fraction  of  excited 

atoms  relaxing  via  the  1th  type  of  radiationless  transition,  hv^  Is  the  enere*. 
of  the  phonons  emitted  In  the  1th  process,  and  is  the  phase  angle  at  which 
the  signal  due  to  transitions  from  the  Initial  state  of  the  1th  transition  Is 
maximum.  Np  and  1,^^  g^e  the  number  and  energy  of  the  photons  in  the  excitation 
beam.  For  the  case  of  Nd3+  Ions  the  sum  In  Eq.  (2)  should  Include  three 
different  types  of  radiationless  decay  processes  which  occur  during  the  relax- 
ation of  the  Ion  from  the  level  excited  by  the  incident  light.  First,  phonons 
will  be  emitted  during  the  decay  from  the  absorption  state  to  the  ^^3/2  ®®t®8table 
level.  Second,  the  ^^2/2  does  not  have  unit  radiative  quantum  efficiency 

and  thus  some  amount  of  relaxation  from  this  state  to  the  ground  state  occurs 
radlatlonlessly.  Third,  the  majority  of  the  fluorescence  transitions  from 
the  ^^2/2  ^^8te  terminate  on  levels  of  the  ^Ij  multlplets  other  than  the 
ground  state  and  subsequent  relaxation  to  the  ground  state  occurs  radlatlonlessly. 

Previous  measurements  have  shown  that  there  Is  virtually  unit  probability  of 
relaxing  to  the  metastable  level  after  pumping  Into  higher  excited  states. 

The  radiative  quantum  efficiency  of  the  ^^3/2  ^8  of  tf*®  order  of  “0.9  In 

the  absence  of  quenching  interactions. ^3  Most  of  the  emission  terminates  on 
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either  the  *•***  ^^9/2  ““itlplets.  Thus  from  these  cons Iderst ions 

end  the  energy  level  diagram  shotm  in  Fig.  2,  it  would  seem  that  the  higher 

lying  absorption  bands  should  be  more  intense  in  the  PAS  spectrum  compared  to 

the  lower  energy  bands  than  in  the  absorption  and  excitation  spectra  due  to 

the  first  type  of  decay  processes  listed  above.  This  is  the  type  of  behavior 

observed  in  other  samples. Figure  10  shows  that  this  is  not  the  case  here 

as  exemplified  by  the  fact  that  the  8000  A band  is  equal  or  greater  In  intensity 
€ 

than  the  5800  A band. 

This  comparison  of  the  absorption,  excitation,  and  photoacoustic  spectra 

indicates  that  there  is  relatively  more  radiationless  quenching  occurring  after 

o • 

excitation  into  the  8000  A band  than  after  excitation  into  the  5800  A band. 

Since  the  absorption  coefficients  of  these  two  bands  are  different,  a different 

spatial  distribution  of  excited  Ions  will  occur  after  excitation  into  each  band. 

For  such  heavily  concentrated  systems  as  those  considered  here,  the  high  absorption 

coefficient  of  the  8000  A band  will  result  in  a significant  fraction  of  the 

excited  Nd  ions  being  located  close  to  the  surface  of  the  crystal.  The  somewhat 

e 

smaller  absorption  coefficient  for  the  5800  A band  allows  for  more  of  the  excited 
Ions  to  be  located  farther  into  the  bulk  of  the  crystal.  One  mechanism  to  ex- 
plain the  spectral  properties  reported  here  Is  surface  quenching  which  is  dis- 
cussed in  detail  in  the  next  section. 

Vr.  INTERPRETATION  OF  RESIH.TS 

At  this  point  it  is  worthwhile  to  summarize  the  experimental  observations 
which  must  be  explained  by  any  viable  model  of  concentration  quenching  in 
NdxYi^_]jP50]^^  and  similar  materials.  The  quenching  rate  must  vary  linearly  with 
Nd  concentration.  Also  it  must  have  a very  weak  dependence  on  temperature 
and  it  must  be  proportional  to  the  absorption  cross  section.  The  quenching 
mechanism  cannot  lead  to  variations  in  oscillator  strengths  with  Nd  con- 
centration, cannot  result  in  spectral  energy  transfer,  and  cannot  lead  to 
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non-exponential  fluorescence  decays.  These  requirements  along  with  quanti- 
tative theoretical  predictions  place  very  rigid  restrictions  on  developing 
an  acceptable  model  for  concentration  quenching  in  these  materials.  The 
mechanism  responsible  for  concentration  quenching  of  Nd  fluorescence  In  other 
types  of  materials  such  as  Y2Al50j^2  •Nd^'*'  is  cross-relaxation  between  pairs  of 
closely  spaced  ions.  This  leads  to  a quadratic  dependence  of  the  quenching 
rate  on  Nd  concentration  which  Is  contrary  to  the  linear  dependence  observed 
for  the  pentaphosphate  hosts.  The  other  common  mechanism  for  concentration 
quenching  which  is  observed  in  various  materials  is  migration  of  the  energy 
among  the  active  Ions  until  it  reaches  a sink  where  It  is  transferred  to  an 
ion  '/hlch  can  dissipate  the  energy  radiatlonlessly.  This  mechanism  is  not 
consistent  with  several  of  the  requirements  mentioned  above.  There  is  no 
dependence  of  the  quenching  rate  on  absorption  cross  section  In  this  model  and 
generally  quenching  by  this  mechanism  has  been  found  to  lead  to  a strong 
temperature  dependence.  A quantitative  estimate  of  the  concentration  of  sinks 
necessary  to  give  the  observed  quenching  rate  can  be  obtained  In  the  fast 
diffusion  limit  from  the  expression  Wq-AirDRgCg  where  D is  the  diffusion  coef- 
ficient. Rfi  Is  the  exciton  trapping  radius  at  the  sink,  and  Cg  is  the  sinx 
concentration.  For  the  case  of  neodymium  pentaphosphate  a sink  concentration 
of  only  10  ppm  of  the  Nd  Ions  is  necessary  to  cause  the  observed  quenching 
rate  whereas  site-selection  spectroscopy  shows  that  significant  numbers  of  Nd 
ions  exist  In  sites  of  different  local  perturbations  without  resulting  in 
energy  transfer  or  quenching. 

Since  the  normal  models  for  concentration  quenching  do  not  appear  to 
apply  to  stochlometrlc  laser  materials  It  Is  necessary  to  develop  a new  model. 

One  model  which  appears  to  meet  all  of  the  requirements  outlined  above  Is 
based  orlmarllv  on  surface  quenching  of  the  excitation  energy.  In  this  model  the 
Nd  ions  on  the  surface  of  the  sample  act  as  sinks  where  the  energy  is  dissipated 
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radlationlessly.  This  gives  a high  concentration  of  sink  sites  but  only  In  a 


very  small  region  and  not  distributed  uniformly  throughout  the  sample.  This 
situation  will  result  in  effective  quenching  for  concentrated  materials  where 
the  excitation  energy  is  deposited  near  to  the  sample  surface  whereas  it  will 
be  much  less  effective  in  lightly  doped  materials  where  the  energy  is  deposited 
more  uniformly  throughout  the  sample.  Similarly  for  a given  sample  quenching 
will  be  more  effective  where  the  absorption  of  exciting  light  is  stronger  due 
to  the  fact  that  more  of  the  excited  ions  will  reside  closer  to  the  surface. 

Weaker  absorption  bands  will  allow  excitation  energy  to  be  deposited  further 
into  the  crystal  resulting  in  less  effective  quenching. 

To  determine  whether  or  not  a surf>:ce  quenching  model  satisfies  all  of 
the  necessary  requirements  for  explaining  the  observed  data,  an  expression  for 
the  quenching  rate  can  be  derived  from  simple  geometric  arguments  based  on  the 
situation  shown  in  Fig.  11.  The  incident  light  excites  a Nd  ions  at  a distance 
R from  the  surface  of  the  crystal  and  the  energy  can  as  an  exclton  migrate  out 
from  this  site  to  other  Nd  ions.  According  to  standard  diffusion  theory, 
after  a time  equal  to  the  intrinsic  lifetime  of  the  exciton  the  excitation 
energy  will  be  found  on  one  of  the  ions  on  the  surface  of  a sphere  of  radius  I 
centered  at  the  position  of  the  originally  excited  ion.  £ is  the  diffusion 
length  which  is  determined  by 

where  Tq  is  the  intrinsic  lifetime  of  the  exclton.  If  it  is  assumed  that  only 
excitons  which  migrate  far  enough  to  reach  the  surface  are  quenched,  the  quenching 
rate  for  ,ar...eKt.ited  ion  at  depth  R is 


Wq  (R)  ■ px”^  (A'/A) 

Wq  (R)  - 0 

where  A is  tne  area  of  the  sphere  of  radius  £,  A*  is  the  part  of  the  area  of 
this  sphere  which  is  outside  the  sample  surface,  and  p is  a factor  to  account 


(4) 
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for  the  probability  of  energy  being  dissipated  radlatlonlessly  Instead  of 
radlatlvely  when  It  reaches  the  surface.  The  fraction  of  the  area  outside  the 
sample  can  be  found  from  the  double  Integration 


2lf  ^dip2lj 


,it/2 


sin  0dO 


m 


- 41  Vji2Ir2  Cos'^(R/t). 

Thus  t\»p  part  of  Eq.  (4)  becomes 

Wq  (R)  - p(TyTr£)-l  y£2.g2  co8‘^(R/£)  0<  R££ 

The  observed  fluorescence  Intensity  Is  proportional  to  the  concentration 
of  excited  Ions.  At  time  t the  remaining  concentration  of  excited  Ions  which 
were  originally  created  at  a depth  R Into  the  sample  Is  given  by  the  product 
of  the  probability  of  exciting  an  Ion  at  that  position  multiplied  by  the  ex- 
ponential decay  function  for, such  Ions 

OrCi)  - exp^-t(T^“^+«p(R^lQoCexp-oCR  dR 

Here  C is  the  concentration  of  Nd  Ions  and  o Is  the  absorption  cross  section 
and  Iq  Is  the  Intensity  of  the  Incident  light.  The  total  concentration  of 
excited  Nd  ions  Is  found  by  Integrating  this  expression  over  the  entire  thick- 
ness of  the  crystal  T.  Using  Eqs.  (4)  and  (6)  this  becomes 
n (t)  - oCIge  [-tp(Tcn)"^  Vi-(R-<-)^  cos~^^  / > J dR 

' + poCRj  dR» 

Evaluation  of  this  expression  shows  that  In  this  model  the  time  dependence  of 
the  ''’tiorescjencc  .int.t.’nsi  Hi  <N‘scrlbed  by 

I(t)  - 1(0)  exp^-t(T^^  + 0.226po  C«,^Q)j 
Thus  the  quenching  rate  Is  given  by 


0. 226pnC"T 


■1 


"Q  • 0 

The  results  of  these  considerations  show  that  a surface  quenching  model 
predicts  exponential  fluorescence  decays  and  a quenching  rate  which  varies 
linearly  with  both  Nd  concentration  and  absorption  cross  section.  Since  the 


(5) 


(6) 


(7) 


(8) 
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ion- Ion  interaction  responsible  for  the  migration  is  a resonant  interaction,  no 
thermal  activation  is  necessary  and  only  a weak  temperature  dependence  is  expected. 
The  quenching  rate  should  actually  be  sometrhat  greater  at  low  temperatures  trtiere 
there  is  less  phonon  scattering  to  limit  the  Man  free  path  of  the  exclton.  This 
appears  to  be  verified  bv  the  lifetime  data  shown  in  Fig.  5.  Quantitatively  the 
experimentally  observed  quenching  rate  for  neodymium  pentaphosphate  can  be  pre- 
dicted by  Eq.  (10)  if  the  thermal  conversion  probability  for  excltons  reaching 
the  surface  (p)  is  taken  to  be  unity  and  the  diffusion  length  is  of  the  order  of 
3,600  A.  Previous  estimates  of  the  diffusion  length  of  excitons  in  this  material 
have  ranged  from  360  to  5,000  A.  A theoretical  estimate  for  the  diffusion  co- 
efficient can  be  obtained  from  the  expression 


D - 1 

6 ” 

A rough  estimate  of  the  hopping  rgte  between  two  Nd.loos  seoprated  bv, a 
distance  can  be  found  by  assuming  a nearest  neighbor  electric  dipole- 
dipole  Interaction  for  two  ions  in  exact  resonance 

^ H ” '^Nd-Nd  ” 

where  the  critical  interaction  •llstan'-e  Ro.i®  given  by 


r 3e^  f^  n T„  "jii 
0 “ |^8tt2  n^  vj,  • 


Here  f is  the  oscillator  strength  of  the  transition  between  the  ground  and 
metastable  state.  V is  the  average  wavenumber  in  the  region  of  spectral  over- 
lap, and  a is  the  spectral  overlap  Integral.  For  two  overlapping  Lorentslan 
lines 


(11) 

(12) 

(13) 


” “ if 

where  v°  is  the  peak  position  and  is  the  Lorentzian  contribution  to  the  line 
width.  Thus  for  perfect  overlap  Q is  proportional  to  the  inverse  of  the  homo- 
geneous contribution  to  the  spectral  line  width.  This  is  measured  at  low 
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Ceaperaturea  to  be  about  6 cm~^.  This  can  be  used  In  Eq.  (13)  which  leads  to 

e 

a predicted  value  of  54  A for  Rq.  Then  for  a nearest  neighbor  distance  of 

a 

5.6  A Eq.  (12)  predicts  an  Ion-ion  interaction  rate  of  3.65  x 109  sec“l.  Sub-  _ 

stituting  these  numbers  into  Eq.  (11)  yields  a predicted  diffusion  coefficient 

of  1.9  X 10~^  cm^  sec~^.  Then  substituting  this  result  into  Eq.  (3)  gives 

o 

a predicted  diffusion  length  of  about  5,000  A which  is  consistent  with  the 
value  needed  to  Interpret  the  concentration  quenching  data  with  a surface 
quenching  model . 

An  attempt  was  also  made  to  fit  the  data  with  a surface  quenching  model 
in  which  the  excitation  energy  did  not  diffuse  but  rather  was  transferred  to 
the  surface  ions  in  a single  step  electric  dipole-dipole  process.  The  pre- 
dictions of  this  model  were  not  consistent  with  the  observed  results  because 

they  predicted  a non-exponential  decay  time  and  quantitatively  an  unphysically  i 

large  value  for  the  interaction  strength  had  to  be  assumed  to  give  the  correct 
magnitude  for  the  quenching  rate. 

A direct  check  was  made  on  the  prediction  that  the  quenching  rate  is 
proportional  to  the  absorption  cross  section  by  monitoring  the  fluorescence 
lifetime  as  the  laser  excitation  was  scanned  from  the  peak  to  the  wings  of 
the  absorption  band.  The  results  obtained  at  14  K for  the  NdP30]^4  sample  are 
shown  in  Fig.  12.  There  is  an  obvious  increase  in  fluorescence  lifetime  and 
thus  a decrease  in  the  quenching  rate  as  the  excitation  is  scanned  from  the 
region  of  high  absorption  to  one  of  small  absorption  as  predicted.  However  the 
dependence  of  Wq  on  o appears  to  be  less  than  the  linear  dependence  predicted 
by  Eq.  (10).  There  may  be  several  reasons  for  this.  First  It  is  difficult 
to  change  the  laser  wavelength  without  also  altering  its  power  and  beam  position. 

Thus  it  is  Impossible  to  exactly  repeat  the  experimental  conditions  from 
point  to  point  in  Fig.  12  and  only  the  general  trend  of  the  data  can  be 
considered  important.  Also  the  theoretical  derivation  of  Eq.  (10)  is  based 
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! I on  the  assumption  that  the  number  of  photons  absorbed  from  the  Incident  beam 

; ! 

at  depth  R is  directly  proportional  to  the  absorption  cross  section.  This  Is 
only  true  for  small  values  of  the  absorption  coefficient.  The  exact  expression 
Is 


n ■ oq  [l-exp(-oCR)3 

1, 

- HqOCR  [ 1->soCR+6(oCR)2-.  . .] 


(15) 


and  thus  for  the  high  optical  densities  treated  here  the  variation  with 
absorption  cross  section  should  be  less  than  linear.  The  linear  approximation 
was  made  In  deriving  Eq.  (10)  only  to  simplify  the  calculations. 

There  is  still  a question  as  to  the  exact  mechanism  for  quenching  after 
the  exclton  has  migrated  to  the  surface.  The  energy  levels  of  the  surface 
Ions  may  be  perturbed  In  such  a way  that  they  act  as  sinks  or  allow  for  cross- 
relaxation to  occur  with  other  Nd  Ions.  A third  possibility  Is  that  there  Is 
such  a high  density  of  excited  Nd  Ions  on  the  surface  that  bl-excltonlc  pro- 
cesses are  responsible  for  the  quenching.  The  latter  possibility  was  Investigated 

by  monitoring  the  fluorescence  lifetime  as  a function  of  laser  power  for  the 

o 

NdP^Oj^^  sample  at  14  K.  The  excitation  was  near  the  peak  of  the  5800  A band 
and  the  power  was  decreased  by  a set  of  calibrated  neutral  density  filters. 

The  fluorescence  lifetime  was  found  to  Increase  slightly  with  decreased  laser 
Intensity.  This  suggests  that  some  amount  of  exclton-exclton  interaction  may 
be  taking  place.  However  the  dependence  of  Wq  on  laser  power  does  not  appear 
to  be  strong  enough  to  account  for  all  of  the  quenching  and  most  probably 
several  different  mechanisms  are  taking  place  at  the  same  time. 


VII.  DISCUSSION  AND  CONCLUSIONS 

The  surface  quenching  model  described  In  the  preceeding  section  Is  con- 
sistent with  all  of  the  results  obtained  In  this  Investigation.  It  also  appears 
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to  be  consistent  with  data  reported  by  other  workers.  Similar  surface  quenching 
effects  have  also  been  observed  In  other  types  of  materials  with  high  absorption 
coefficients  such  as  aromatic  hydrocarbon  crystals. 

For  crystals  with  high  absorption  coefficients  radiative  reabsorption 
effects  should  be  expected.  We  attempted  to  determine  the  extent  of  such 
effects  by  observing  the  fluorescence  lifetime  as  a function  of  sample  align- 
ment. The  lifetime  appeared  to  lengthen  slightly  In  going  from  front  face  to  back 
face  Illumination.  All  data  reported  here  were  obtained  with  front  face 
Illumination  to  minimize  radiative  reabsorption  effects. 

It  should  be  pointed  out  that  there  Is  a possible  alternative  explanation 
for  the  fact  that  the  relaxation  processes  In  the  excited  states  do  not  appear 
to  contribute  to  the  PAS  signal.  This  would  be  true  If  the  decay  occurred 
radlatlvely  instead  of  radiatlonlessly.  There  has  been  some  weak  fluorescence 
observed  between  1.4  and  1.7  pm  which  has  not  been  explained.^  However  this 
was  observed  only  at  a high  level  of  excitation  and  only  at  room  temperature 
and  thus  can  not  generally  account  for  all  of  the  observed  quenching  character- 
istics. 

It  Is  still  not  clear  exactly  why  the  fluorescence  quenching  rate  Is  so 
much  less  In  NdP^O]^^  than  In  YAG:Nd.  One  obvious  difference  Is  that  there  is 
less  resonance  between  the  ^^3/2  ^^15/2  transitions  and  the  ^19/2  * 

^^15/2  ^^Botptlon.  Also  a distinguishing  feature  of  the  pentaphosphate  host  is 
that  the  Nd  Ions  are  fairly  isolated  from  each  other  and  do  not  share  the  same 
oxygen  ion.^  Another  Interesting  observation  Is  that  the  quenching  of  Nd  flu- 
orescence In  phosphate  glass  hosts  is  also  significantly  less  than  for  other 
glasses. Which  of  these  features  Is  most  Important  In  determining  the  quenching 
characteristics  Is  not  evident.  It  Is  known  that  hydrogen  and  other  Impurity  ions 
are  very  effective  in  quenching  the  fluorescence  In  these  materials. These 
Impurities  may  perturb  the  energy  levels  of  neighboring  Nd  ions  to  allow  stronger 
quenching  Interactions  to  occur.  The  energy  levels  of  surface  Ions  may  similarly 
be  perturbed. 
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interactions  to  occur.  The  energy  levels  of  surface  ions  may  similarly  be 
perturbed. 

The  numerical  factors  in  Eq.  (3)  and  (11)  are  exact  for  Isotropic  random 
walks  on  a cubic  lattice.  These  will  be  slightly  different  for  the  penta- 
phosphate  lattice  structure  but  the  cubic  lattice  approximation  is  sufficient 
for  the  order  of  magnitude  estimates  needed  here. 

The  diffusion  coefficient  of  the  order  of  ''<10“^  cm^  sec”^  determined  for 
NdPjOj^^  is  significantly  greater  than  the  typical  values  of  D which  range  be- 
tween 10“^  to  10“^^  cm2  sec”^  for  energy  migration  among  other  systems  of  rare 
earth  ions  in  solids. 27  The  reason  for  this  is  primarily  the  higher  concentration 
of  Nd  ions.  This  results  in  efficient  spatial  diffusion  without  spectral 
diffusion.  As  the  excitation  energy  migrates  near  to  a spectrally  different 
site  there  are  always  a significant  number  of  ions  nearby  in  equivalent  sites 
that  favor  transfer.  Only  at  the  surface  is  there  a significant  concentration 
of  perturbed  ions  to  cause  effective  quenching. 

In  summary  the  results  of  time-resolved  site  selection  spectroscopy, 
photoacoustic  spectroscopy,  and  other  investigations  presented  here  are  con- 
sistent with  a model  involving  spatial  energy  migration  without  spectral  dif- 
fusion which  results  in  significant  radiationless  quenching  at  the  surface. 
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A)  Block  diagram  of  time-resolved  spectroscopy  setup 

B)  Block  diagram  of  photoacoustic  spectroscopy  setup 

Energy  level  scheme  of  Nd^"*" 

Excitation  Spectrum  of  ^F'3/2  level  In  NdQ  cYq  sPeO^^  and 
NdPjOi^  at  300  K.  • • 3 

Concentration  dependence  of  the  ^F3/2  fluorescence  lifetime  and 
quenching  rate  for  NdjjY2^^jjP50j^^  at  14  K 

Temperature  dependence  of  the  fluorescence  lifetime  of  the^^Fj^2 
level  in  NdP50. . for  laser  excitation  wavelengths  of  5804  A 
and  5836  %. 

^^3/2  '*  ^^9/2  fluorescence  in  NdP^O^^  20  ps  after  selective 
excitation  into  different  regions  of  the  ^G,.y»  absorption 

band  at  14  K 1 01 ^ 

Line  1 of  ^^3/2  ^^9/2  fluorescence  in  NdQ  ^Yg  5P50j^^  after 

20  us  excitation  at  various  wavelengths. 

A)  Time  dependence  of  the  fluorescence  deviation  function  after 
excitation  at  5784.5  X and  5837  X in  NdP50j^^  for  temperatures 
14,  42,  75  and  110  K 

B)  Time  dependence  of  thg  fluorescence  deviation  function  after 
excitation  at  5805.5  A and  5835  A in  NdQ.sYQ  5P50j^^  for 
temperatures  14  and  77  K 

Temperature  dependence  of  t^e  fluorescence  duration  function 
after  excitation  at  5784.5  A and  5837  X in  NdP5034 

Photoacoustic  Spectra  of  Ndo,9Yo.iP5034  and  NdPjOj^^  at  300  K 

Top  and  side  view  geometry  for  excitation  energy  migration  in 
surface  quenching  model. 

^^3/2  fluorescence  lifetime  vs.  wavelength  of  laser  excitation 
at  14  K along  with  laser  probed  absorption  band  at  300  K for 
NdP50i4 
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